Modelos Multinivel

Modelos mixtos y de partición de varianza

Plot example

#------------------------------------------------
# power analysis function
#------------------------------------------------

get_power <- function(n_size, d){
library(dplyr)
power <- pwr::pwr.t.test(n=n_size,d=d,sig.level=.05,alternative="two.sided") %>%
         purrr::pluck('power')

table <- data.frame(n = n_size, power = power, d = d)

return(table)
}


#------------------------------------------------
# generate data
#------------------------------------------------

table_n <- dplyr::bind_rows(
purrr::map_df(as.list(seq(from = 10, to = 500, by = 1)),~ get_power(n_size = ., d = .1), .id = 'index'),
purrr::map_df(as.list(seq(from = 10, to = 500, by = 1)),~ get_power(n_size = ., d = .2), .id = 'index'),
purrr::map_df(as.list(seq(from = 10, to = 500, by = 1)),~ get_power(n_size = ., d = .3), .id = 'index'),
purrr::map_df(as.list(seq(from = 10, to = 500, by = 1)),~ get_power(n_size = ., d = .4), .id = 'index'),
purrr::map_df(as.list(seq(from = 10, to = 500, by = 1)),~ get_power(n_size = ., d = .5), .id = 'index'),
purrr::map_df(as.list(seq(from = 10, to = 500, by = 1)),~ get_power(n_size = ., d = .6), .id = 'index'),
purrr::map_df(as.list(seq(from = 10, to = 500, by = 1)),~ get_power(n_size = ., d = .7), .id = 'index'),
purrr::map_df(as.list(seq(from = 10, to = 500, by = 1)),~ get_power(n_size = ., d = .8), .id = 'index'),
purrr::map_df(as.list(seq(from = 10, to = 500, by = 1)),~ get_power(n_size = ., d = .9), .id = 'index'),
purrr::map_df(as.list(seq(from = 10, to = 500, by = 1)),~ get_power(n_size = ., d = 1.00), .id = 'index')
)

Attaching package: 'dplyr'

The following objects are masked from 'package:stats':

    filter, lag

The following objects are masked from 'package:base':

    intersect, setdiff, setequal, union

#------------------------------------------------
# display plot
#------------------------------------------------

library(ggplot2)
ggplot(table_n, aes(d, power, group=n, colour=n)) + 
geom_line() +
geom_hline(yintercept = .80, lintype = 'dashed', color = 'red') +
theme_minimal() +
theme(
legend.position = 'right',
# legend.title = element_blank(),
text = element_text(size = 12, family = 'Avenir'),
panel.background = element_blank(),
panel.grid.minor.y = element_blank(),
panel.grid.major.y = element_blank(),
panel.grid.minor.x = element_blank(),
panel.grid.major.x = element_blank(),
axis.text.x =  element_text(colour = "grey30", size = 10),
axis.text.y =  element_text(colour = "grey30", size = 10),
axis.line = element_line(colour = "black")
) +
xlab('effect size') +
ylab('satistical power')

Warning in geom_hline(yintercept = 0.8, lintype = "dashed", color = "red"):
Ignoring unknown parameters: `lintype`

Test titles

A second subtitle

And A third subtitle